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I. INTRODUCTION 

The experimental analysis of EEG waveforms has 
been dominated by spectral decomposition and ampli- 
tude analysis, mainly autocorrelogram techniques 
though crosscorrelation technique has been used for 
phase measurements. These techniques as applied to 
the analysis of EEG signals are abundant in litera- 
ture. To cite a few of them, we note that the ear- 
liest ones were done by Brazier and Barlow [1], [2], 
[3], [4], and notable simplifications done by Kamp, 
et al . [5], DeBoer and Kuyper [6], and Vo-Ngo4 et 
al [7], The fact that the phases of EEG signals as 
compared to a reference clock signal contain infor- 
mation about the mental alertness has been pointed 
out to us by Anliker [8] who has used phase-vector 
and contour-graphic techniques to study the phase 
entrainment phenomenon. Adey and Walter [9] have 
used a phase detection technique in the analysis of 
EEG records in the cat. The EEG a-rhythm has been 
treated as the result of mutual synchronization of 
a population of spontaneously oscillatory processes 
by Wiener [10]. The potential usefulness of the 
phase entrainment in either recognizing the patho- 
logical states of the brain or classifying mental 
states might be widely explored if investigators had 
a reliably simple direct measurement technique. In 
this paper we describe a frequency discrimination 
technique as realized by an inexpensive frequency 
discriminator for direct on-line measurement of the 
entrainment phenomenon in EEG as entrained by the 
frequency of a sensory stimulus. We demonstrate the 
use of this device for detecting the presence or ab- 
sence of the stimulus effect and the measurement of 
the time delays in the entrainment. 

In our study, the EEG signal is first filtered 
by a narrow-band filter with center frequency about 
the a-rhythm of the individual. We shall denote 
this filtered EEG signal as alpha signal a(t). It 
is then reasonable to consider the alpha signal as 
narrow-band random process with the alpha frequency 
as the mean frequency of the spectral band. Then 
a sample function of this random process is express- 
ed as [11] 

a(t) = a (t) cos [2Tr%t+*(t)] (1) 

If the bandwidth of its power density spectrum is 
much smaller than its mean frequency Iq, then the 
processes a(t) and 't’(t) in Equation (1) will be 
slowly varying functions of time as compared to 
cos 2TiiQt so that the interpretation of a(t) and 
<t>(t) as envelope and phase has meaning. The alpha 
signal process a(t) can thus be interpreted as 
amplitude and phase modulated signal process with 


carrier frequency fj-j . The process 4(t) can also 
be interpreted as a relative phase angle, in the 
sense that the alpha signal process a(t) differs 
in phase from the signal (cos by <t> { t) . It 

is only in this manner that we are able to ascribe 
meaning to the phase of EEG as a time-varying quan- 
tity. The fact that the angular velocity or fre- 
quency in radians per second is the time derivative 
of the angular position leads us to define the in- 
stantaneous frequency ff(t) in cycles per second 
or Hertz by 


f. (t) 

l 


1 d <t> 
2 tt dt 


( 2 ) 


The second term in Equation (2) which is proportion- 
al to the rate of change of the time-varying phase 
<t>(t) can be interpreted as the instantaneous fre- 
quency deviation relative to the a frequency IJ-j 
provided it is stable. Our technique to be describ- 
ed in detail provides a direct measurement of (d*)/ 
(dt) . The study on the statistical properties of 
the random processes a(t) and *(t) will be re- 
ported in another article. 

When the brain is excited by an outside stimu- 
lus frequency f s and the EEG is entrained by the 
stimulus, then we should have 

f (t) = f = f t + T < t < t + T (3) 
i Cts nn —i x 


where t n and t f signify respectively the time 
instants at which the stimulus is on and the stimu- 
lus is off, and T n denotes the time delay between 
t n and the time instant at which the entrainment 
occurs and Tf, the time delay between tf and the 
time instant at which the entrainment disappears. 
Equation (3) implies that 


d$ 

dT 


= 0 , 


t +T < t < t„ + T- 
n n — — f f 


(4) 


Since our technique measures (d<t>)/(dt) directly, 
the measurement of EEG signals which were recorded 
while the human subjects with eyes closed were stim- 
ulated by stroboscopic flashes at the rate of 10 
flashes per second for one minute, then no stimula- 
tion for one minute, then another minute of stimula- 
tion, etc., does show that (dil>)/(dt) = 0 for al- 
ternate one minute intervals. 

The schemes used for the measurement of (d*)/ 
(dt) , T n , and Tf , and for the automatic detection 
of the states of stimulus-on and stimulus-off in the 
EEG are to be described next. 


II. METHOD 

An ideal frequency discriminator should pro- 
duce an output voltage linearly dependent on input 
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frequency. There are many ways for the realization 
of frequency discrimination. The most commonly used 
one is the so-called balanced demodulator with the 
well-known S curve. The scheme which we have used 
here is the zero-crossing detector type as described 
in [12]. In order to eliminate the influence of the 
amplitude variation a(t), we use an amplifier fol- 
lowed by a hard limiter in front of the detector. 

The output of the frequency discriminator is propor- 
tional to the instantaneous frequency f^ (t) as ex- 
pressed in Equation (2), and does not depend on the 
amplitude of the signal. We are, however, interest- 
ed in the rate of change of <t>(t) . The value fj^ 
in Equation (2) is subtracted. In our device fj-j 
can be preset at any value in the range of 5-15 Hz. 
The device thus yields (K/2n) (db/dt) on-line. The 
characteristic curve of the device is given in Fig.l. 
In Fig. 2, we show the functional block diagram of 
the analysis scheme. Detailed circuitry of the de- 
vice is given in Fig. 3. 

The response time or delay T n (or tf) is 
defined as the length of time elapsed since the on- 
set (or end) of stimulus to the time instant at 
which the brain wave is entrained (or desynchroniz- 
ed) . This delay T n (or Tf) is, in a sense, the 
time' which the brain takes to synchronize (or desyn- 
chronize) . There is a noticeable delay between the 
onset of stimulus-on (or stimulus-off) and the time 
at which the entrainment is measurable. This fact 
is shown in Fig. 4 where the (d40/(dt) as the out- 
put of the frequency discriminator does not reach a 
flat plateau (or start to wander) immediately after 
the onset of stimulus-on (or stimulus-off) . The 
fluctuation of the x(t) ^ (db)/(dt) curve is mea- 
sured by the quantity 

P = !> I X <V - ra J <5) 

i=l 

within a specified window width N, where mu is 
the mean value of (d<t>)/(dt) or x(t) in the 
window; i.e., 

N 

m = i y x(t.) (6) 

N N l 

1=1 

New values of p are calculated by sliding the win- 
dow. The onset of the flatness (or the wandering) 
is determined by the time instant at which p is 
below (or above) a preset threshold that is deter- 
mined experimentally. The response time or delay 
T n (or Tj) is measured as the time lapse between 
the onset of stimulus-on (or stimulus-off) and the 
onset of the flatness (or the wandering) of (db)/ 

(dt) as measured by p. Owing to the delay intro- 
duced by the window width, a correction constant is 
subtracted from the above obtained time intervals. 

The detection scheme utilizes the same program 
as that used in measuring the response time or de- 
lay, but without the prior knowledge of the onset 
of stimulus-on and stimulus-off. In other words, 
the influence of stimulus (entrainment) is said to 
be detected whenever the quantity p reaches a val- 
ue below a preset threshold . 

These schemes are realized by digital computer 
programs in PAL III language used on a PDP-8 compu- 
ter. Both of these measurements can be performed 
on-line. Figure 5 depicts the flow chart of the 
computer program for the response-time measurement. 

III. RESULTS 

Although the measurements were performed on EEG 
signals recorded on magnetic tape, they can be per- 


formed on-line. As mentioned before, these EEG sig- 
nals were recorded while the subjects with eyes 
closed were stimulated by stroboscopic flashes for 
one minute at the rate of 10 flashes per second, 
then no stimulation for one minute, then another min- 
ute of stimulation, etc. The flashes were generated 
by Grass photo stimulator model PS-2E on its inten- 
sity scale No. 2. The EEG signals are the differen- 
tial potentials between the electrodes placed at 
left parietal and left occipital. The EEG signals 
are referred to as either stimulus-entrained or non- 
entrained . In order that the EEG signal be more ap- 
propriately described by the narrow-band-process mod- 
el, the EEG signal was passed through a narrow-band 
filter with a band-width of 1 . 5 Hz centered at 10 Hz 
to obtain the alpha signal a(t) as described pre- 
viously . 

A section of the output of the frequency dis- 
criminator along with the same section of EEG sig- 
nals of subject B being analyzed and a square wave 
indicating when the stimulus was turned on or off 
are shown in Figs. 6 and 7. The display of apparent 
high frequency components in the raw EEG of Figs. 6 
and 7 are due to the slower writing speed (30 cm per 
minute) of our strip chart recorder as compared to 
the conventional writing speed. It is clearly seen 
that there is a delay following the onset of the 
stimulus before entrainment is detected; the EEG re- 
mains synchronized for a short period of time also 
after the stimulus was turned off. In Figs. 8, 9, 
and 10, we show the output of the frequency discrim- 
inator with the early, middle, and late sections of 
EEG signals of subject H as the input along with a 
square wave indicating the events of stimulus-on and 
stimulus-off. The offset of time coordinates in the 
above figures is due to the positions of pens of the 
strip chart recorder. From these figures, we observe 
that subject H synchronized with the stimulus very 
well at the beginning section (13th min. to 19th 
min.) of the experiment as shown by the measurement 
of (d<t)/(dt) depicted in Fig. 8. As shown in Fig. 
9, however, the subject's EEG signal did not entrain 
to the stimulus frequency as well during the middle 
(105th min. to 111th min.) of an experimental ses- 
sion which lasted about three hours. There is no 
sign of entrainment of the EEG signal in the end 
section (176th min. to 182nd min.) of the experiment 
as shown in Fig. 10. The capability to synchronize 
with the stimulus may be a measure of certain brain 
states such as alertness. Further study is needed 
for this measurement by this technique. The fre- 
quency discrimination technique as realized by an 
inexpensive device is effective in measuring the en- 
trainment phenomenon of EEG signals. 

It is possible that variations in the delays 
T n and Tf in synchronization and desynchroniza- 
tion, respectively, might give some indication of 
change of levels of alertness, although our data 
does not include alertness estimates. Results of 
the measurement of delays for both subjects are 
shown in Figs. 11 and 12, respectively. The abscis- 
sas in the figures denote the time since the begin- 
ning of the experiments and the ordinates are the 
measured delays. Any T n > 25 seconds is considered 
as not synchronized. Subject B seems to have longer 
response time to stimulus-on in the first hour of 
experiment than those in the second hour. Subject H 
did not entrain so well after the first hour of ex- 
periment. In general, we see that it takes longer 
for the subjects to synchronize with the stimuli 
than it takes to return to the natural states after 
the stimuli were shut off. The statistics of t n 
and tf and their correlation with the mental 
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states of subjects need further study. However, the 
technique for measuring "these delays is satisfactory. 

In Figure 13, we show a typical section of the 
conditions or states of the stimulus-on and stimu- 
lus-off of EEG as results of computer analysis of 
the output of the frequency discriminator in com- 
parison with the timing track for the on and off of 
the stimulus. This was generated by a double- 
threshold detector. The thresholds were determined 
experimentally. We analyzed the records of both 
subject B and subject H. The ratio of the sum of 
the differences between the time of the on and off 
conditions of EEG as determined by the detector and 
on-off timing track for the stimulus to the total 
time of on and off as measured by the timing track 
was calculated. This ratio is 28.6 per cent for 
subject B and 35.8 per cent for subject H. We did 
not take into account the time delays introduced by 
the smoothing window and the frequency discrimina- 
tor in the evaluation of these ratios. The values 
of these ratios may be reduced. The use of the 
above-mentioned ratio as a performance measure of 
the detector is appropriate if we do not have the 
knowledge of the length of time and the manner in 
which the stimulus is turned on or off. With the 
knowledge that the stimulus was turned on for one 
minute, then shut off for one minute, then was on 
for another minute, etc., and counting only the 
number of the on-off cycles as determined by the 
detector, we obtain the following result: 

(1) For subject B, there are 57 on-off cycles detec- 
ted by the detector as compared to 57 on-off cycles 
on the timing track (a 100 per cent detection) . 

(2) For subject H, there are 56 on-off cycles detec- 
ted by the detector as compared to 59 on-off cycles 
on the timing track (a 94.7 per cent detection). 

IV. CONCLUSIONS 

The frequency discrimination technique describ- 
ed here has been effective in measuring the entrain- 
ment phenomenon in EEG signals. We have shown that 
the narrow-band-process model used here for charac- 
terizing signals is valid. The technique as real- 
ized by an inexpensive device described previously 
provides an effective way for further study on the 
synchronization phenomenon as related to alertness. 
Statistical properties of the phase ♦(t) must be 
gathered for further interpretation of mental states. 
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FIG. 1 FREQUENCY DISCRIMINATOR OUTPUT CHARACTERISTIC. 


FREQUENCY DISCRIMINATOR 
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FIG. 2 BLOCK DIAGRAM OF THE ANALYSIS SCHEME. 



FIG. 3 FREQUENCY DISCRIMINATOR - CIRCUIT DIAGRAM. 



FIG. 4 RESPONSE TIME OR DELAY MEASURED. 



FIG. 5 FLOW CHART FOR 
THE RESPONSE-TIME MEASUREMENT PROGRAM. 
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EEG INPUT TO THE SCHEMATIC DIAGRAM SHOWN IN FIG. 2 
WHEN STIMULUS IS ON. STRIP-CHART RECORDER WRITING 
SPEED IS 30 CM PER MINUTE. 



EEG INPUT TO THE SCHEMATIC DIAGRAM SHOWN IN FIG. 2 
WHEN STIMULUS IS OFF. STRIP-CHART RECORDER WRITING 
SPEED IS 30 CM PER MINUTE. 



FIG. 8 MEASUREMENT OF -J- OF THE ALPHA SIGNAL OF THE 
at 

EARLY PART OF EEG RECORD OF SUBJECT H . 



ON OFF 


FIG. 9 MEASUREMENT OF ~ OF THE ALPHA SIGNAL OF THE 
dt 

MIDDLE PART OF EEG RECORD OF SUBJECT H. 



OFF ON 


FIG. 10 MEASUREMENT OF — OF THE ALPHA SIGNAL OF 
dt 

THE LATE PART OF EEG RECORD OF SUBJECT H. 


DETECT ON 
OUTPUT 



FIG. 13 THE STATES OF STIMULUS-ON AND STIMULUS-OFF 
OF EEG AS RESULTS OF THE ANALYSIS IN COMPARISON 
WITH THE ACTUAL ON AND OFF OF THE STIMULUS • 
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Fig. 11. RESPONSE TIME STUDY OF SUBJECT B. 




Fig. 12. RESPONSE TIME STUDY OF SUBJECT H. 
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